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A STEEPEST DESCENT METHOD 
FOR OSCILLATORY RIEMANN-HILBERT PROBLEMS 



P. Deift and X. Zhou 



In this announcement we present a general and new approach to analyzing the 
asymptotics of oscillatory Riemann-Hilbert problems. Such problems arise, in par- 
ticular, in evaluating the long-time behavior of nonlinear wave equations solvable 
by the inverse scattering method. We will restrict ourselves here exclusively to the 
modified Korteweg de Vries (MKdV) equation, 

y t - 6y 2 y x + Vxxx = 0, -oo < x < oo, t>0, 
y(x,t = 0) = y (x), 

but it will be clear immediately to the reader with some experience in the field, 
that the method extends naturally and easily to the general class of wave equations 
solvable by the inverse scattering method, such as the KdV, nonlinear Schrodinger 
(NLS), and Boussinesq equations, etc., and also to "integrable" ordinary differential 
equations such as the Painleve transcendents. 

As described, for example, in [IN] or [BC], the inverse scattering method for 
the MKdV equation leads to a Riemann-Hilbert factorization problem for a 2 x 2 
matrix valued function m = m(-; x, t) analytic in C\R, 

/•js m + (z) = m-{z)v Xtt , zeR, 

m(z) — > / as z — > oo, 

where 

m±(z) = limm(z ± ie; x, t), 

eJ.0 
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If yo{x) is in Schwartz space, then so is r(z) and 

r(z) = —r(—z), sup \r(z)\ < 1. 

From the inverse point of view, given v(z), one considers a singular integral equation 
(see [BC]) for the associated quantity fi(z: x, t) — m + (z; x, t)(b+ l ) x . t = m^(z; x, t)(bZ l ) x ,t 
and the solution of the inverse problem is then given by 




where 

w x . t = {w+) x ,t + (w-) x ,t, w± = ±(b± - I). 

Significant work on the long-time behavior of nonlinear wave equations solvable 
by the inverse scattering method, was first carried out by Manakov [M] and by 
Ablowitz and Newell [AN] in 1973. The decisive step was taken in 1976 when 
Zakharov and Manakov [ZM] were able to write down precise formulae, depend- 
ing explicitly on the initial data, for the leading asymptotics for the KdV, NLS, 
and sine-Gordon equations, in the physically interesting region x — 0(t). A com- 
plete description of the leading asymptotics of the solution of the Cauchy problem, 
with connection formulae between different asymptotic regions, was presented by 
Ablowitz and Segur [AS], but without precise information on the phase. The as- 
ymptotic formulae of Zakharov and Manakov were rigorously justified and extended 
to all orders by Buslaev and Sukhanov [BS 1-2] in the case of the KdV equation, 
and by Novokshenov [N] in the case of NLS. 

The method of Zakharov and Manakov, pursued rigorously in [BS] and in [N], 
involves an ansatz for the asymptotic form of the solution and utilizes techniques 
that are somewhat removed from the classical framework of Riemann-Hilbert prob- 
lems. In 1981, Its [I] returned to a method first proposed in 1973 by Manakov in 
[M] , which was tied more closely to standard methods for the inverse problem. In 
[I] the Riemann-Hilbert problem was conjugated, up to small errors which decay 
as t — > oo, by an appropriate parametrix, to a simpler Riemann-Hilbert problem, 
which in turn was solved explicitly by techniques from the theory of isomonodromic 
deformations. This technique provides a viable, and in principle, rigorous approach 
to the question of long-time asymptotics for a wide class of nonlinear wave equa- 
tions (see [IN]). Finally we note that in [B], Buslaev derived asymptotic formulae 
for the KdV equation from an exact determinant formula for the solution of the 
inverse problem. 

In our approach we consider the Riemann-Hilbert problem (1) directly, and by 
deforming contours in the spirit of the classical method of steepest descent, we show 
how to extract the leading asymptotics of the MKdV equation. In particular for 
x < 0, let ±zq = zty / \x\/12t be the stationary phase points for i(4iz 3 + xz). Then 
the first step in our method is to show that (1) can be deformed to a Riemann- 
Hilbert problem on a contour X of shape (see Figure 1), in such a way that the 
jump matrices v x . t onRcS and on the compact part 

of S\R away from ±zo, converge rapidly to the identity as t — > oo. Thus we are 
left with a Riemann-Hilbert problem on a pair of crosses T, A U X s (see Figure 2). 
As t — > oo, the interaction between E A and X s goes to zero to higher order and 
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Figure 1 



Figure 2 



Figure 3 

the contribution to y(x,t) through (2) is simply the sum of the contributions from 
Tj A and S B separately. Under the scalings z — > z(48tzo) -1 ^ 2 =F zo, the problems on 
T, A and S s are reduced to problems on a fixed cross, with jump matrices which 
are independent of time, and which can be solved explicitly in terms of parabolic 
cylinder functions, as in [I]. Substitution in (2) yields, finally, the asymptotics for 
y(x,t). 

Our result is the following: let 



7T 1 f Z " 

<j>(z ) = argT(» — — — argr(zo) + - / log|s - z |d(log(l - |r(s)| 2 ) 
4 71 J-z a 



(here T is the standard gamma function) and let 

/ \ 1/2 

f» = Ut- cos(16tzl - v log(192te3) + 4>(zo)), 



\ 3tzo t 

where v = -(2T:)- 1 log(l - |r(z )| 2 ) > 0. Set r = tz$ = (\x\/12t 1 ^ 2 . 

Theorem. Let yo{x) lie in Schwartz space with reflection coefficient r(z). As 
t — > oo, the solution y(x,t) of MKdV with initial data yo(x), has uniform lead- 
ing asymptotics conveniently described at fixed t » 1, in t/ie six regions shown in 
Figure 3. 

In region I, for any j, 

y(x, t) = y a + 0((-x)- j + (-xy^Cji-x/t)) 

where Cj(-) is rapidly decreasing. In region II, 

y(x,t) = y a + (tz )- 1 / 2 O(r- 1 / 4 ). 
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In region III, 

y(x, t) = (3i)- 1A Mz/(3t) 1/3 ) + 0(r 2 / 3 /t 2/3 ), 

where p is a Painleve function of type II. 
In region IV, 

y(x, t) = (3i)- 1/3 ^/(3t) 1/3 ) + 0(r 2 / 3 ). 
In region V, for any j, 

y (x,t) = (my^pix/iM) 1 / 3 ) + o(r j + r 2 ^ e - 12 ^ 2/3 ) 

for some r\ > 0. 

Finally, in region VI, for any j, 

y(x,t) = 0((x + t)- j ). □ 

Remark 1 . The reader may check that in the overlap regions the asymptotic forms 
do indeed match. Also the reader may check that in regions II and IV, the formulae 
for the leading asymptotics agree with those in [IN] . 

Remark 2. The above error estimates are not the best possible and in region II in 
particular, the t -1 / 4 decay can certainly be improved. 

Remark 3. There is no obstacle in the method to obtaining an asymptotic expansion 
for y{x, t) to all orders. 

Remark 4. As noted at the beginning of this announcement, the method we have 
presented extends naturally to the general class of nonlinear wave equations solv- 
able by the inverse scattering method. Also, there is no difficulty in incorporating 
solutions with solitons. 
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